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PROCEDE DE DETERMINATION D'UNE STRUCTURE EN 3b A PARTIR 
D^UNE mA&E EN 2D, NOTAMMENT LA STRUCTURE D*UN OS 

La ppesente invention concerne un procede de determination d'une structure en 
3D a partir d'une image en 2D, notamment celle d'un os. 

L'osteoporose est une maladie qui affecte en grand nombre surtout les femmes 
apres la menopause vers Tage de 50 ans sachant que cette maladie peut off ecter 
5 toute personne et a tout age. Elle est caracterisee par une faible masse osseuse 
et une deterioration du tissu osseux. Ces degradations conduisent a des risques 
importants de fractures particulierement de la hanche, de la colonne vertebrale 
et du poignet. :j 
II existe bien sur un cevtQm nombre de f acteurs de risques mais ceci ne suff it 
10 pas pour determiner la probabilite pour une personne de se f racturer uh element 
de son squelette sachant qu'il n'y a aucun symptome prealable. 
Cette maladie est courante dans la population et devrait aller en s'etendant du 
fait de la sedentarisation notamment et du vieillissement de la population. 
Or le traitement des personnes off ectees apres fracture est tres long et tres 
15 couteux car ce sont des soins a long terme. Les consequences sont souvent 
graves car conduisant a des invalldites voire a la mort des patients, 
Aussi est-il necessaire de prevenir ce type d'affection mais a condition de 
determiner le risque de fracture, ce qui est du seul ressort du praticien. 
Si Ton arrive a predire les risques, le proticien dispose alors de traitements par 
20 hormonotherapie et par prise d' out res produits de la pharmacopee comme la 
calcitonin ou les biphosphonates, en plus de conseils d' hygiene de vie pour tenter 
d'enrayer cette affection. 
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II existe des appareillages pour determiner la densite osseuse appeles 
densitometres. 

On peut penser que des lors que Ton connoTt la densite osseuse, en comparant 
avec des courbes etalon, ie praticien peut alors determiner les risques et etablir 
un diagnostic adapte et prescrire un traitement satlsf aisant. 
C'est ce qui est couramment utilise. 

Par la demande de brevet US 5 774 520, on sait que la probabilite d'une 
fracture est directement liee a la densite osseuse. 

Les densitometres disponibles dans Ie commerce utilisent les rayons X. ou les 
rayons gamma de faible energie. L'absorption etant proportionnelle d la densite 
de Tos, on peut done penser d correler les deux pour avoir un resultat 
sat isf aisant permettant des comparaisons. Ce brevet prevoit un travail avec 
deux sources d' emission de photons en sorte de traiter la dispersion engendree 
par Ie fait que les matieres osseuses peuvent etre considerees comme 
homogenes mais que les autres matieres liquides, les tissus musculaires, les 
cartilages ne peuvent etre considerees comme une matiere unique, homogene. 
Un autre brevet US 6 385 283 utilise la densite mais combine cette mesure avec 
des images permettant de determiner les risques de fracture. Ces images sont 
realisees sur la colonne vertebrale du patient et Toperateur determine la 
presence d'amorces de fractures, be plus, comme dans Tart anterieur connu, 
ces informations sont comparees avec des donnees acquises au fur et d mesure 
sur de tres nombreux cas. Ce diagnostic est off ine par Ie praticien d partir des 
antecedents du patient et par introduction des facteurs de risques 
supplementaires. 

La demande de brevet WO 86/07531 propose de realiser une image d'un os 
donne, par exemple Ie calcaneum, d un age auquel Ie patient presente 
normalement sa plenitude de masse osseuse, Ensuite, Ie precede consiste d 
realiser des images successives de ce meme os pour Ie praticien de realiser des 
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comparaisons et de determiner I 'evolution pour conclure sur le degre de risque 
de fracture. 

Toutes ces methodes et autres dispositif s sent bases sur la seule mesure de la 
densite de la mat i ere osseuse. 
5 Or pour que le praticien arrive a un diagnostic satisfaisant, il faut qu'il ait des 
informations autres que la seule densite. 

La definition de Tosteoporose est donnee par i* Organisation Mondiale de la 
Sante : desordre caracterise par une faible masse osseuse et des alterations 
de la microarchitecture du tissu osseux, conduisant a augmenter la fragilite de 

10 tos et en consequence a augmenter le risque de fracture " 

C'est ainsi qu'il existe des analyses de matiere osseuse in vivo plus sophistiquees 
a partir d'analyses de sang et d'urines faisant intervenir des marqueurs osseux 
mais ces analyses sont couteuses, contraignantes pour le patient et de toutes les 
foQons donnent des indications de composition de la matiere mais pas sur ,sa 

15 Structure. 

Or, on peut comprendre que la matiere osseuse presente une resistance 
differente suivant qu'il y a plus ou moins de matiere en quantite ce qui resulte 
de la mesure de densite. Par contre, pour une meme densite osseuse done pour 
une meme quantite de matiere, suivant T architecture de la structure osseuse, la 

20 resistance mecanique peut varier grandement, 

Ce qui est utile pour un praticien, ce sont les proprietes mecantques de Tos d'un 
patient induisant directement les risques de fracture. Quand il y a diminution de 
la resistance osseuse, il y a conjugaison de deux f acteurs, la baisse de la densite 
minerale osseuse DMO et la desoiteration de la structure osseuse. 

25 A ce jour, il n'existe pas de moyens permettant de determiner la structure 
osseuse. On peut utiliser de Timagerie par resonance magnetique ou le scanner 
mais ce sont des examens lourds et couteux, surtout s'il faut effectuer un suivi 
longitudinal des patients. 
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De plus, la capture d' images en 3b a parti'r de moyens d'imagerle est pour 
{'instant limitee par les capacites meme de ces nr\oyens car la resolution sjjatiale 
est superieure aux dimensions de la plupart des travees de liaison permettant de 
realiser la connectivite, c'est-a-dire les liaisons entre les noeuds du reseau en 
5 trois dimensions. 

Le prpcede selon la presente invention est maintenant decrit en detail, afin de 
permettre par I'analyse d'une image 2D d'estimer les parametres mecaniques. 
Selon I'invention, le procede de determination du coefficient de resistance 
ultime Cu de la partie trabeculaire d'un os a partir d'une image 2D numerisee, 

10 obtenue par Imagerie utillsant un rayonnement X ou un rayonnement gamma de 
faible energie, se caracterise essentiellement en ce que Ton realise une 
correlation entre la densite minerale osseuse DAAO determinee d partir de cette 
image 2D et un parametre de structure obtenu a partir de la meme image 2D. 
De fagon preferentielle, le parametre de structure retenu est la connectivite 

15 osseuse Xv- 

Selon une autre caracteristique relative au meilleur mode de realisation connu au 
jour du depot, le procede consiste a determiner un parametre a lie a la 
connectivite osseuse Xv obtenu par la succession des etapes suivantes : 

- choisir au hasard un pixel de I 'image 2D qui est a un niveau de gris h(0), 

20 - choisir une droite initiale ayant une direction determinee egalement au hasard, 

- se deplacer d'une distance r et determiner le niveau de gris h(r) de ce nouveau 
point, 

- calculer la variance des niveaux de gris avec la formule : V(r) = [h(r) - h(0)]^, 

- tracer la courbe associee log-log, et 

25 - determiner la pente de cette courbe log-log qui represente ledit parametre a, 
Le procede consiste a realiser une correlation entre la densite minerale osseuse, 
DMO, obtenue a partir de I'image 2D et ledit parametre a en sorte de 
determiner la contrainte ultime Cu de la partie trabeculaire de Tos concerne. 
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Selon unc autre revendication, la correlation preferentielle retenue dans le cas 
de la presente invention est du type exponentielle ; 
Cu = bo + bi* exp (bz * DMO) a 

- avec bo, bi, bz qui sont des coefficients obtenus par regression non 

5 lineal re. 

Dans la presente invention, on assimile la densite minerale osseuse bMO a la 
fraction volumique de soiide fs= / Vo dans laquelle : 

- Vo : volume d'analyse (mm^), et 

- Vs > volume du soiide (mm^) 

10 et le parametre a d la densite de connectivite Xv = X / Vo, en mm^^ 
avec la relation x = po + Pe - Nep 

- po * nombre d'amas de la partie soiide, ce nombre est generalement 

egal d 1. 

- P2 • nombre de surfaces internes, generalement egal d 0. 
15 - Nep • nombre d^Euler-Poincarre 

Pour supporter cette description, des dessins sont annexes et les figures qu'ils 
presentent, illustrent la description et montrent essentiellement des resultats 
obtenus. 

- figure 1 : courbe des resultats avec un modele lineaire, en f onction de Cu' 

20 - figure 2 : courbe des resultats avec un modele exponentiel, Cu en f onction de 

- figure 3 : representation graphique de Cu en f onction de Xv 

- figure 4 : representation de f s en f onction de Xv 

- figure 5 :courbe de la f onction V(r) 

25 - figure 6 : determination du parametre a d partir de la courbe log-log de la 
f onction V(r) 

- figure 7 : courbe bMO en f onction de f s 

- figure 8 : courbe de a en f onction de Xv 
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- figure 9 t courbe DMO en f onction de a 

- figure 10 : modele lineaire Cy en f onction de DMO 

- figure 11 : courbe de en f onction 

L*etude porte sur un os trabeculaire ou spongieux qui constitue environ 207o de 
5 la matiere osseuse. Vos cortical qui Tentoure assure la rigidite de Tos complet 
et done du squelette tandis que I'os trabeculaire assure la reprise et la 
resistance aux forces de con\pression, 

Selon r invention, on considere que la structure trabeculaire comprend des 
travees qui sont interconnectees et la resistance mecanique de T ensemble 
10 resulte de cette connectivite, c'est-a-dire du nombre de boucles fermees. La 
resistance mecanique resulte aussi du degre de mineralisation de ces travees, 
mineralisation qui est un parametre important. 

Le procede consiste a analyser une image numerisee obtenue a partir d'une 
emission de rayons X. Cette image est une projection sur une surface d'une 
15 Structure 3D ayant ete traversee par les photons emis. 

Chaque pixel qui constitue T image peut etre analyse independamment mais alors 
on ne peut determiner que la densite minerale osseuse DMO si bien qu'il convient 
non seulement d' analyser un pixel de f agon isolee mais de les analyser les uns par 
rapport aux autres. 

20 Le spectre du rayonnement doit etre le plus monochromatique possible afin 

d'evlter une trop forte derive dans Tenergie des photons. 

En effet, si Tenergie varie la penetration variera et donnera des nuances de gris 
sur r image pouvant induire des erreurs en faisant croire a une absorption la ou il 
n*y en a pas. 

25 A partir de cette image obtenue avec un spectre etroit, on sait que chaque 
volume elementaire subit une action directe des photons, perpendiculairement 

mais aussi une action indirecte des rayons incidents diffuses par les volumes 
elementafres voisins apres qu'ils aient ete eux-memes ivawerses. 
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Af in de supprimer la part de diffusion, on peut utiliser de preference deux 
energies distinctes. On dispose de 2 equations a 2 inconnues et on peut done 
deduire la part d' attenuation reelle qui nous interesse dtrecten\ent et la part 
resultant de la diffusion, 
5 Si Ton recherche la contrainte ultime Cu (MPa) , il faut tenir compte des 
parametres suivants ; 

- Vo • volume d' analyse (mm^) 

- Vs volume du solide (mm^) 

- fs » fraction volumique de solide 
10 avec la relation fs = Vs / Vo 

- Po : nombre d amas de la partle solide (c'est-d-dire Tensemble des parties 
connectees ou disconnectees), ce nombre est generalement egol a 1. 

- p2 • nombre de surfaces internes, c'est-d-dire les trous resultant de la micro 
porosite Merne des travees or a Techelle de travail, cette porosite ne se vqit 

5 pas, ce nombre est generalement egai d 0. 

- Nep : nombre d' Euler-Poincarre 

- X * connectivite 

avec la relation x = Po + Pa - Nep 

- Xv densite de connectivite (mm"^) 
20 avec la relation Xv = X / Vo 

Le procede selon la presente invention vise a permettre de correler deux 
parametres mesurables et de trouver la relation qui les lie. 
1 / Si Ton utilise un modele lineaire tel que : 
25 Cq' = ao'^fs-^-ai 

- Cu' etant une valeur predite de Cu 

- oo et Qi sent des coefficients de regression lineaire, 
il s*agit Id d'une simple correlation de Cuovec fs. 
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La courbe obtenue est representee sur la figure 1 et Von constate que la 
dispersion est tres importante. Si Ton se place par exemple a fs constant, 
Tecort des valeurs de Cu est important. 

On ne peut retenir ce modele lineaire dont les resultats sont insuffisants en 
5 precision mais qui peut neanmoins servir de comparatif , 

2/ Si I 'on utilise un modele exponentiel tel que ; 

Cu* = bo + bi* exp(b2 ^ fs) Xv 
Avec bo, bi, hz qui sont des coefficients obtenus par regression non lineaire 
10 (methode des moindres carres generalisee). 

La dispersion est f ortement limitee, voir figure 2. est estime a partir de fs 

Ces variations sont representees sur la figure 3. 

Si fset Xv augmentent, Cu augmente. Cecl est parfaitement naturel, si la fraction 
15 volumique de solide augmente et la connectivite augmente, alors la contrainte 
ultime augmente. 

Si on fixe Xv alors Cu est directement fonction de fs, de fagon exponentielle. A 
connectivite constante, la contrainte ultime augmente en fonction de 
r augmentation de la densite osseuse. 
20 Si on fixe fs alors Cy est directement fonction de xv. de fagon lineaire mais la 
contrainte ultime Cudiminue lorsque la connectivite augmente. 
Si la fraction solide est constante, on constate que la resistance ultime Cu 
diminue lineairement en fonction de la connectivite. 

Plus on a de ncEuds dans une structure et plus cette structure devient fragile, 
25 ceci de fa^on surprenante, contrairement a un prejuge bien etabli. 

On peut alors tracer le modele de representation de la figure 4 qui montre les 
valeurs de Cu. 
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Pour la determination de la connectivite, 11 faut done trouver un parametre qui 
suive une meme loi et qui soit le plus independent possible de la densite osseuse 
pour ne pas etre influence. 

On pourrait realiser un calcul a partir de deux parametres mesurables comme la 
5 fraction volumique osseuse et la connectivite, ces deux parametres etant 
obtenus notamment par Imagerie a Resonance Magnetique, 

Par contre, Tobtention des valeurs par cette methode met en oeuvre des 
appareillages couteux, peu disponibles, ce qui conduit a des predictions d"un prix 
de revient eleve interdisant un suivi longitudinal et done une diffusion massive 

10 aupres des patients. 

II faut done pouvoir trouver deux valeurs mesurables par exemple a partir d'une 
simple densitometries qui est, elle, tout a fait abordable et qui peut etre 
renouvelee regulierement. En effet, des densitometries sont deja realisees de 
fagon massive pour de nombreuses applications. La diffusion aupres des patients 

15 peut etre tres important e. 

Le precede selon la presente invention consiste d determiner ces deux 
parametres mesurables comme la fraction volumique osseuse et la connectivite, 
ou plus exQctement des parametres qui pourraient leur etre correles de fa^on 
suf f isamment etroite pour donner des resultats satisfaisants. 

20 Le premier parametre rapidement mesurable et fiable est la densite minerale 
osseuse, Des appareils ont ete developpes pour eff ectuer ces mesures de fagon 
fiable et reproductible. Si Ton trace la courbe DMO en fonction de fs. on 
constate que la dispersion est faible et que Ton peut utiliser DMO en lieu et 
place de fs,. Voir figure 7. Le premier parametre est obtenu. 

25 Le second parametre est appele a et determine a partir d'une image 2D 
numerisee par exemple obtenue avec un rayonnement X. 

Pour determiner ce nouveau parametre a, on choisit tout d'abord au hasard un 
pixel de T image 2D qui est a un niveau de gris h(0). 
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On choisit ensuite une droite isiitide ayant une direction determinee egalement 
au hasard. 

On se deplace d'une distance r et on determine le niveau de gris h(r) de ce 
nouveau point, 

5 On calcule alors la variance des niveaux de gris, ce qui donne : 

Le processus est iteratif et aleatoire pour determiner a et permet d'obtenir 
par calcul un ensemble de points et la courbe associee est representee sur la 
figure 5. On fait ainsi converger la statistique V(r). 
10 Ainsi, plus on determine de points ayant un meme niveau de gris a proximite du 
pixel choisi, plus on se trouve dans la matiere et done plus on est homogene et 
plus on doit trouver de ncBuds de connexion. 

Si Ton trace la courbe log-log de cette fonction, figure 6, on trouve sur les cinq 
premiers points en T occurrence, un comportement lineaire et la pente de cette 
15 droite est le coefficient denomme a et propre a la presente invention. 

Si Ton trace maintenant la courbe a en fonction de la connectivite volumique Xv. 
figure 8, on constate qu'il y a une correlation suffisante pour que Ton puisse 
rempiacer Xv par a. 

On peut deduire le diagramme de la figure 9 de ces courbes precedent es, 
20 diagramme qui montre la relation entre DMO, a et Cu , 

DMO et a sont deux parametres directement mesurables et/ou pour le moins 
calculables a portir d'une meme image en 2D. 

On constate sur ce diagramme que lorsque la denstte minerale osseuse reste 
constante, la contrainte ultime diminue lorsque a augmente. 
25 Ceci est bien equivalent au modele exponentiel precedemment indique avec la 
fraction osseuse en fonction de la connectivite, voir figure 4, sauf que la 
connectivite Xv n'etait pas mesurable directement alors que a est calcule a 
partir de r image 2b, done quantifiable. 



1 er depot 



11 

On peut done ecrire la meme relation nriais avec Ies parametres determinables a 
partir de la meme image : 

Cu = bo + bi'^ exp(b2 * DMO) * a 
C'est la representation de la figure 9. la courbe de la figure 11 caracterise ce 

5 modele exponentiel de Cu' en fonction de a. On constate une faible dispersion 
compare au modele lineaire dont la courbe est mentionnee h la figure 10, 
Ceci est un exemple non limitatif qui montre la possibilite de lier la densite 
osseuse dont la valeur est determinee d partir d'une image 2b, a un parametre 
de structure determine a partir de cette meme image 2b pour obtenir la 

10 contrainte ultime Cu. 
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REVENlblCATTONS 

1. Precede de determination du coefficient de resistance ultime Cu de la partie 
trabeculaire d'un os a partir d'une image 2b numerisee, obtenue par Imagerie 
utiiisant un rayonnement X ou un rayonnement gamma de faible energie, 
caracterise en ce que Ton realise une correlation entre la densite minerale 

5 osseuse DMO determinee d partir de cette image 2D par tout moyen adapte a ce 
type d'image et un parametre de structure obtenu a partir de la meme image 2D. 

2. Procede de determination du coefficient de resistance ultime Cu selon la 
revendication 1, caracterise en ce que le parametre de structure retenu est la 
connectivite osseuse Xv. 

10 3. Procede de determination du coefficient de resistance ultime Cy selon la 
revendication 2, caracterise en ce qu'il consiste a determiner un parametre a lie 
d la connectivite osseuse xv obtenu par la succession des etapes suivantes : 

- choisir au hasard un pixel de 1' image 2D qui est a un niveau de gris h(0), 

- choisir une droite initiale ayant une direction determinee egalement au hasard, 
15 - se deplacer d'une distance r et determiner le niveau de gris h(r) de ce nouveau 

point, 

- calculer la \/ariance des niveaux de gris avec la formule ; V(r) = [h(r) - h(0)3^, 

- tracer la courbe associee log-log, et 

- determiner la pente de cette courbe log-log qui represente ledit parametre a. 
20 4. Procede de determination du coefficient de resistance ultime Cu selon la 

revendication 3, caracterise en ce que Ton realise une correlation entre la 
densite minerale osseuse, DMO, obtenue a partir de cette image 2D et ledit 
parametre a en sorte de determiner la contrainte ultime Cu de la partie 
trabeculaire de Tos concerne. 
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5. Procede de determination du coefficient de resistance ultime Cy selon la 
revendication 4, caracterise en ce que la correlation est du type exponentielle : 

Cu' = bo bi^^ exp (b2 ^ DMO) ^ a 

- avec bo, bi, be qui sont des coefficients obtenus par regression non 

5 lineaire. 

6, Procede de determination du coefficient de resistance ultime Cu selon la 
revendication 4 ou 5, caracterise en ce que l*on ossimile la densite minerale 
osseuse DMO d la fraction volumique de solide fs= Vs / Vq dans laquelle : 

- Vo • volume d' analyse (mm^), et 
10 - Vs : volume du solide (mm^) 

et le parametre a d la densite de connectivite Xv = X / Vo, en mm"^ 
avec la relation x ^ Po + Pz - Nep 

- po • nombre d'amas solides, ce nombre est generalement egal d 1. , 

- p2 ' nombre de surfaces internes, generalement egal d 0. 
15 - Nep • nombre d'Euler-Poincarre 
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